Association of molecules using a resonantly modulated magnetic field 
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We study the process of associating molecules from atomic gases using a magnetic field modulation that is 
resonant with the molecular binding energy. We show that maximal conversion is obtained by optimising the 
amplitude and frequency of the modulation for the particular temperature and density of the gas. For small 
modulation amplitudes, resonant coupling of an unbound atom pair to a molecule occurs at a modulation fre- 
quency corresponding to the sum of the molecular binding energy and the relative kinetic energy of the atom 
pair. An atom pair with an off-resonant energy has a probability of association which oscillates with a fre- 
quency and time-varying amplitude which are primarily dependent on its detuning. Increasing the amplitude 
of the modulation tends to result in less energetic atom pairs being resonantly coupled to the molecular state, 
and also alters the dynamics of the transfer from continuum states with off-resonant energies. This leads to 
maxima and minima in the total conversion from the gas as a function of the modulation amplitude. Increasing 
the temperature of the gas leads to an increase in the modulation frequency providing the best fit to the thermal 
distribution, and weakens the resonant frequency dependence of the conversion. Mean-field effects can alter the 
optimal modulation frequency and lead to the excitation of higher modes. Our simulations predict that resonant 
association can be effective for binding energies of order h x 1 MHz. 

PACS numbers: 03.75.Nt, 34.50.-s, 34.20.Cf 
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I. INTRODUCTION 



Cold diatomic molecules are often produced from atomic 
gases by varying a magnetic field around a zero-energy reso- 
nance lUsyt] or by photoassociation Jj]. Generally applica- 
ble implementations of the former technique are linear ramps 
of the magnetic field across the resonance HI d, 0, 0, H, Isl Q3. 
muni J3l [Till , and fast switches to fields close to the reso- 
nance [ 1 ]. In addition, the long lifetime of cold 6 Li 2 dimers al- 
lows them to be created by holding the magnetic field close to 
resonance, causing thermalisation of the atomic into a molec- 
ular gas fl5l|. In recent experiments, dimers of 85 Rb2 lfl6ll and 
85 Rb- 87 Rb [17] were associated from cold atomic gases by ap- 
plying a magnetic field modulation resonant with the molec- 
ular binding energy. This technique eliminates the need for 
the magnetic field to spend time in the near-resonant, strongly 
interacting region. It therefore reduces the unwanted effect 
of heating llq, fl8ll during the production of molecules. The 
narrow Fourier spectrum of the pulse accurately targets the 
molecular state, minimising the coupling to deeper bound 
states and highly energetic continuum states. In a direct com- 
parison to a linear ramp using the same apparatus, Thompson 
et al. reported more efficient conversion using resonant asso- 
ciation tig] . This technique has also been used as an accu- 
rate probe of molecular binding energy Il6lll7ll . In addition, 
radio-frequency pulses have been used to associate 40 K- 87 Rb 
dimers IU9I1 . 

In this paper we study the resonant association of molecules 
from thermal and condensed gases. Our approach precisely 
accounts for the continuum of states in a gas. The transi- 
tion amplitude from a pair of unbound atoms to the bound 
molecular state depends on the relative kinetic energy of the 
atom pair. A resonant continuum energy exists, at which the 
transition amplitude to the molecular state increases linearly 
with time. At small modulation amplitudes, the resonant con- 
tinuum energy is given by the difference between the energy 



corresponding to the modulation frequency and the molecular 
binding energy. The distribution of atoms in different contin- 
uum states, all contributing to the molecular production, gives 
the total conversion a dependence on temperature. The width 
of the thermal distribution increases with the temperature of 
the gas, weakening the resonant behaviour of the molecu- 
lar production. The continuum distinguishes the current case 
from the association of atom pairs held in optical lattices l2(ill . 
where the resonant modulation couples the discrete ground 
state of the tightly confining potential to the molecular bound 
state. 

We find damped oscillations in the number of molecules 
produced in the short-time limit, as observed in Ref. flrill . The 
damping is caused by the dephasing of the transition ampli- 
tudes from states across the continuum. After the damping 
out of the initial oscillations, the conversion increases at a rate 
which displays resonant dependence on the modulation fre- 
quency. Maximal conversion is achieved when the frequency 
and amplitude of the magnetic field modulation are together 
optimised for the temperature and density of the gas. The 
modulation amplitude required depends on the sensitivity of 
the molecular state to the magnetic field. As the modula- 
tion amplitude increases states of lower continuum energy are 
resonantly coupled, until the zero-energy continuum state is 
reached. Beyond this point all continuum states are coupled 
in a non-resonant manner. There remain some modulation am- 
plitudes where, for momenta close to the peak of the thermal 
distribution, the transition amplitude is large enough to lead to 
a revival in conversion efficiency. Our calculations of molecu- 
lar production for binding energies ranging from hx5 kHz to 
h x 1 MHz predict that resonant association can be effective 
over this range. We also examine short pulses in pure conden- 
sates, where the mean-field shift and the excitation of higher 
modes alter the dynamics. In condensates, the damping of os- 
cillations in conversion due to the dephasing of the transition 
amplitudes from different continuum states is suppressed. 
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In Sec. [II] we introduce the magnetic field profile used for 
resonant association, and also set up the notation used in our 
calculations. We then discuss the dynamics of the transition 
amplitude for a pair of atoms to a molecule, and the depen- 
dence this has on the continuum, in Sec. [En] In Sec. [IV] 
we examine in turn the effects of altering the duration, fre- 
quency and amplitude of the modulation on the efficiency of 
the molecular production. We also discuss the dependence 
of the conversion efficiency on the temperature and density 
of the atomic gas. In each section we discuss the results of 
Thompson et al. lfl6tl . which formed the original motivation 
for our studies, and then consider resonant association under 
a broader range of conditions. We conclude in Sec.lVl 




II. MAGNETIC FIELD SEQUENCE 



t (ms) 
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In this section we introduce the magnetic field sequence 
used in resonant association, as well as the notation used in 
this paper. Our calculations of molecular production from 
thermal gases use a two-channel approach J2lll22ll23il24ll25ll . 
In the implementation of Ref. 12611 . the two-channel, two-body 
Hamiltonian for the case of a time-varying magnetic field B( t) 
is given by 

H 2B (B(t)) = |bg)i/ bg (bg| + W|bg)(cl| 

+ |cl)(bg|W + |cl>// cl (B(f))<cl| . (1) 

Here, W is the interchannel coupling between the entrance 
and closed channel spin configurations which are indicated by 
'bg' and 'cl', respectively. Choosing the zero of energy to 
coincide with the dissociation threshold of the entrance chan- 
nel makes the closed-channel Hamiltonian contain all of the 
magnetic field dependence of H2v{B(t)). We make the single- 
resonance approximation 1 26], neglecting all closed-channel 
states which are far detuned from E — 0. The single closed- 
channel state retained is referred to as the resonance level 
|0 res ), and is degenerate with the entrance-channel dissocia- 
tion threshold at the field strength B res , i.e. E KS (B KS ) = 0. 
The closed-channel Hamiltonian is then given by 

W)) = l0 reS >K+£™ d (O]<0re S |. (2) 

Here E£ = ^(B„ - B res ), E™ d (t) = [B(t) - B„], B, v is 

the average magnetic field during the pulse, and -^f is the dif- 
ference in magnetic moment between the entrance and closed 
channels. The measurable location Bq of the singularity in 
the scattering length a is shifted from B re . s by the interchannel 
coupling. 

A schematic of the magnetic field sequence used in the ex- 
periments of Thompson et al. [16] is shown in Fig.Q] which 
also summarises our notation. Gases of 85 Rb atoms in the 
(F = 2,;«f = -2) excited Zeeman state were prepared at 
162 G. Experiments were performed using thermal gases of 
temperature T = 20 - 80 nK, as well as partially and wholly 
condensed gases. Following a 5 ms ramp to an average field in 
the range 156 - 157 G, a sinusoidal magnetic field pulse was 
applied for a duration of up to 38 ms. The pulse resonantly 



FIG. 1: (Color online) A schematic of the magnetic field sequence 
used in Ref. 1 16]. The ramps before and after the resonant pulse were 
shown not to create molecules. We assume a jump in 6(f) at t f when 
B(t f ) * fi av . 



associated molecules in the highest excited vibrational bound 
state, which we term Feshbach molecules. The magnetic field 
was then ramped back to 162 G in 5 ms, and the gas held un- 
til all molecules were lost due to spin relaxation 1114 
Absorption imaging before and after this sequence showed the 
depletion of the atomic gas. 

In our calculations of molecular conversion from thermal 
gases we neglect the 5 ms ramps on either side of the pulse, 
which were found in Ref. not to lead to molecular pro- 
duction. We have verified that including the ramps in our 
calculations causes only a negligible difference to the final 
result. The magnetic field pulse is taken to have the form 
B(t) = B av + B mod sin(w mod f), where fi mod and « mod are the 
amplitude and angular frequency of the modulation, as shown 
in Fig.Q] B(t) is assumed to satisfy B(0) - B(t { ) = B dv . Where 
B(tf) + B dw , a sudden jump at ff returning to B av is assumed, 
as illustrated in Fig.Q] Consequently, we analyze the initial 
and final states in terms of the eigenstates of ff™ = H2b,(Bm). 
The local density approximation is valid for weakl y co nfin- 
ing traps such as that used by Thompson et al. Ill6ll28ll . and 
for simplicity we assume a homogeneous gas. The eigenstate 
describing the continuum state of a free pair of atoms with rel- 
ative momentum p is the dressed scattering state \<p™) satisfy- 
ing H™^)™) - (/? 2 /m)|0p V ), where m is the atomic mass. The 
boundary conditions of |0p V ) correspond to an incoming plane 
wave plus an outgoing spherical wave I26ll29tl . The Feshbach 
molecular state at B av , |^ v ) , satisfies Hf^tfig) = E™\dg), 
where E™ = E\,(B m ) is the molecular bound state energy 112611 . 

III. DYNAMICS OF AN ATOM PAIR 

In this section we discuss the two-body dynamics of a pair 
of atoms in the presence of a continuum, due to the mag- 
netic field pulse introduced above. Numerically integrating 
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the Schrodinger equation associated with the Hamiltonian of 
Eq. (0 for a particular pulse gives the two-body time evolu- 
tion operator f/2Bfe 0), linked to H2s(B(t)) by 

d 

in-U 2B (t, t') = H 2B (B(t))U2B(t, f) . (3) 
at 

After the pulse, the wavefunction of a pair of atoms that were 
initially in the state \<p™) has an overlap with the Feshbach 
molecular state given by the transition amplitude 

T(p,t f ) = (<t>l v \U 2 B(.t 1 ,0)\<?>v)- (4) 

This in turn gives the probability density for a transition be- 
tween a free pair of atoms of relative momentum p and a Fesh- 
bach molecule to be 

p(p,t { ) = \T(p,t { )\ 2 ■ (5) 

In this paper we consider resonances where a spherically sym- 
metric resonance level is coupled to the entrance channel by 
spin exchange (3J . Consequently, the transition amplitude and 
probability density depend only on the modulus of the mo- 
mentum. The transition probability density gives the dynam- 
ics of only one state in the continuum. This is to be distin- 
guished from the conversion itself, which includes the contri- 
butions of all the continuum states. 

In the limit of short times and a small modulation ampli- 
tude, U2s{tf, 0) may be approximated using time-dependent 
perturbation theory. Treating the oscillating component of the 
Hamiltonian of Eq. (0 as a perturbation to H'^ gives the first 
order approximation to the two-body evolution operator: 

+ i j ' dt U™(tf - f)l0 le s> J E™° d (O<0 1 -esl C/|b(0 • 

(6) 

Here, the two-body evolution operator at B av is given by 
t/jJgC?) = exp(-iH^t/K). Projecting the estimate of Eq. © 
onto (0^ v | on the left and |0p V ) on the right gives an approxi- 
mation to the transition amplitude of Eq. 

T m (P^ = ~^B mod e-<'^C(p) 
2n oB 

x \ e uo + t t sinfrMf) _ „^ sin(fc>_f f ) l _ (?) 

Here 

u ± = (E™ + fuo mod - p 2 /m) /(2») , (8) 

and C(p) = (<^ v |0 res )(0re S |0p V ) is the product of the overlaps 
of the resonance state with the bound and scattering states at 
B av . Since Z?? v , ftco mo< i and p 2 /m can all be of the same order 
of magnitude, it is not in general possible to make the rotating 
wave approximation and neglect the <y_ term in Eq. 0. 

As implied by the first-order estimate of Eq. 0, the fastest 
growth in the transition probability density p(p, t) for small 
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FIG. 2: (Color online) Evolution of the transition probability density 
p{p, t) in a 85 Rb gas for near-resonant continuum energies. The res- 
onance condition of Eq. (|9j is fulfilled at p 2 /m = h x 0.73 kHz. The 
solid lines show the numerical results, and the dashed lines show the 
perturbative estimate of Eq. (0. The transition probability densities 
shown are labelled with the energy of the continuum state in h x kHz. 
The inset shows the evolution of p(p, i) for the 0.98 kHz contin- 
uum state at longer times. Here, B av = 156.45 G, 5 mo d = 0.065 G, 
Egjh = -5.77 kHz, and a> mo J2n = 6.5 kHz. 



Bmod occurs for the resonant continuum energy p 2 es /m, which 
satisfies 

ZC+fi«W- — =0. (9) 
m 

This corresponds to the sum of the relative kinetic energy 
of the atom pair and the molecular binding energy |£? v | be- 
ing exactly matched by the modulation frequency. As shown 
in Fig. [2] the growth in the transition probability density 
p(p,t) for p 2 /(mh) = 0.73 kHz is quadratic. This corresponds 
to the transition amplitude of Eq. having a linearly increas- 
ing amplitude, similar to a resonantly driven harmonic oscil- 
lator. 

At sufficiently short times, p(p, t) grows quadratically for 
states detuned from the resonant continuum energy. Figure [2] 
shows that p(p, t) displays behaviour different to p{p KS , f) af- 
ter a time of order ft/(\p 2 - p 2 es \/m). The inset of Fig.|2]shows 
the oscillatory nature of p(p, t) , reproduced at short times by 
the analytic estimate of Eq. (0. The numerical approach, ac- 
counting precisely for the continuum of states, yields an enve- 
lope in the oscillation amplitude as well as a frequency shift 
from the estimate of Eq. 0. The thermal gases measured in 
Ref. flql had temperatures in the range 20 - 80 nK. The energy 
ks x 50 nK corresponds to h x 1 kHz. Consequently, the phase 
difference between the continuum states spread over the ther- 
mal distribution becomes significant after times of the order 
of milliseconds. 
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IV. EXPERIMENTAL PARAMETERS AFFECTING THE 
CONVERSION EFFICIENCY 

In this section we study the variation of the molecular pro- 
duction with the duration, frequency and amplitude of the 
magnetic field modulation, and the temperature and density 
of a thermal or fully condensed gas. We refer to the fraction 
of atoms converted to molecules as the conversion efficiency. 
In the limit of small depletion of a thermal atomic gas, this 
is given by a weighted average of the transition probability 
density pip, tf) over a Maxwell distribution: 

— — =2n{2n1\y\ — afpexp \pip,tf) . (10) 

N \mnl J \ m ] 

Here, /3 = l/k^T, n is the density of the atomic gas, N is 
the initial number of atoms and N mo \ is the final number of 
molecules. In this limit the conversion efficiency is propor- 
tional to the density of the atomic gas. We note that this 
approach does not lead to saturation of the conversion effi- 
ciency, which would require the inclusion of genuinely many- 
body effects. This requires the solution of non-Markovian 
Boltzmann-like equations, whose Markov limit has previously 
been used to study the special case of saturation of molecular 
production from magnetic field ramps lf30h . 



A. Pulse duration 

1. Thermal gas 
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FIG. 3: (Color online) Conversion efficiency from a thermal 85 Rb 
gas as a function of pulse duration, for a density of n = 10" cm" 3 . 
All other parameters are the same as those of the data represented 
in Fig. [2] which has been thermally averaged according to Eq. d 1 Ob 
to give the conversion efficiency. The inset shows damped oscilla- 
tions in the conversion efficiency, visible for T = 20 nK but washed 
out for 50 nK by the dephasing of the transition amplitudes from dif- 
ferent continuum states. The dotted lines show the results given by 
thermally averaging the perturbative estimate of Eq. 10. 

efficiency, and so the initial coherence in pip, t) across the dis- 
tribution is destroyed more quickly. 



We first consider the conversion efficiency from a thermal 
gas. The averaging of Eq. (fTTJb gives a contribution from 
p{p, f) for each p, weighted according to the thermal distri- 
bution. Figure [3] shows the resulting conversion efficiency for 
gases of 20, 50 and 80 nK as a function of pulse duration. The 
resonance condition of Eq. (0 is fulfilled at a continuum en- 
ergy of h x 0.73 kHz, which corresponds to 37 nK. Of the gas 
temperatures quoted in Ref. 11611 . 20 nK gives the highest con- 
version efficiency because the most atom pairs have energies 
close to the resonant continuum energy. 

In the experiments of Ref. Ill 611 . damped oscillations in the 
conversion efficiency as a function of time were observed over 
the first few milliseconds. In our calculations, damped oscil- 
lations are visible over approximately 2 ms for a temperature 
of 20 nK. We have verified for several values of B dv and w mo[ j 
that the frequency of the damped oscillations, / con v, for the 
thermal gas case is close to p^JiynK). This is the value of 
a» + /27r for p = in Eq. ([8j, corresponding to the detuning of 
the zero momentum state from the resonant continuum energy. 
Increasing the temperature causes a negative shift in the fre- 
quency of the damped oscillations, together with faster damp- 
ing. For 50 and 80 nK gases our calculations do not predict 
damped oscillations large enough to be observed. The main 
cause of damping is the variation in the oscillation frequency 
of pip, t) with p, as shown in Fig.|2]and discussed in Sec. Hill 
A wider thermal distribution corresponds to a wider spread 
in momentum of the atom pairs contributing to the conversion 



2. Condensed gas 

The critical temperature reported in Ref. [16] is 14 nK, 
with average densities of order 10 11 cm" 3 . For a density of 
10 11 cm" 3 , which we use in our thermal gas calculations, and 
a magnetic field of 156.45 G, the dilute gas parameter Vna 3 
is 0.02. For such gases, which are close to condensation or 
partially condensed, there will be a mean field shift in the fre- 
quency of the oscillations in conversion efficiency. We have 
analyzed this effect for the case of a pure, homogeneous con- 
densate. For our studies of condensed gases we use the cumu- 
lant approach lf3~ll l32l l33ll . In this approach, the atomic mean 
field ^(f) is given by a non-Markovian, nonlinear Schrodinger 
equation Il32ll : 

m-mt) = H m v{i)-y*{i) d T y 2 (T)— h(t,r). (ii) 
at Jo or 

The coupling function hit, r) contains the exact two-body dy- 
namics, and is given by 

h{t, t) = (2;r/0 3 <0| V(f)t/ 2 B(f, T)|O)0(f - t) . (12) 

Here, the evolution operator U2g.it, t) is defined in Eq. (0), |0) 
is the plane wave of zero momentum, V(f) is the diatomic po- 
tential, and Bit - f) is the step function, yielding 1 for t > r 
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and otherwise. The molecular conversion is given by a 
molecular mean field, 

¥ b (r) = -— dT^ 2 (r)— h b (t,r), (13) 
V2 Jo 9t 

where the bound state coupling function is given by 

h b (t,T) = (2^)^ 2 <0 b |f/ 2B (f, T )\O)0(t - t) . (14) 

The coupling functions of Eqs. (TTZl i and (TT~4-b have been de- 
termined from the single-channel approach of Ref . j32tl . The 
magnetic fields used in this calculation are within the range 
for which single-channel approaches have been shown to be 
valid for the 155 G resonance of 85 Rb (H- 
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FIG. 4: (Color online) Frequency of the oscillations in conversion 
efficiency at short times for a pure 85 Rb condensate as a function of 
the initial density. The mean-field shift lowers the frequency of the 
oscillations in conversion efficiency from that given by the two-body 
approach, which is recovered in the limit of low density. The oscil- 
lations are, however, much clearer than those in a thermal gas due 
to the suppression of the contributions of different continuum states 
to the molecular production. The inset shows the variation in time 
of the conversion efficiency for densities of 10 10 (x50 for clarity), 
4 x 10" and 8 x 10" cnT 3 . The 0.5 ms ramp from B = 157.45 G 
to S av = 156.45 G, not shown here, gives a density-dependent initial 
phase to the oscillations. Here B mo *i = 0.065 G, /h = -5.86kHz 
and a) moli /2n = 7 kHz. Because a single-channel approach is used in 
the condensed gas case, the bound state energy is slightly different to 
that given by the two-channel, thermal gas calculations above. 

The oscillation frequency of the conversion efficiency at 
short times, / conv , has a mean-field shift, as shown in Fig. [4] 
In the low-density limit, the value of / conv expected from the 
two-body picture is recovered. The oscillations in conversion 
efficiency are clearer and have weaker damping than those in 
thermal gases, as shown in the inset of Fig. [4] This is due 
to suppression of the dephasing between the transition am- 
plitudes from different continuum states. The main cause of 
damping in this case is the decay of the condensate and molec- 
ular populations into the continuum. In these calculations we 
have included a 0.5 ms ramp from B av + 1 G to B av , in analogy 



to the ramp shown in Fig.Q] Neglecting the ramp and simu- 
lating only the pulse corresponds to instantly turning on the 
interactions at the beginning of the pulse. For the parameters 
used here, this results in strong excitation of higher modes. 
The ramp reduces but does not completely eliminate the ex- 
citations, which are visible in the inset of Fig. [4] as high fre- 
quency oscillations whose amplitude increases with conden- 
sate density. We have extracted the frequency of the damped 
oscillations using the fit procedure of Claussen et al. 13511 . 
which includes exponential damping of the oscillations and a 
linear decay. Strong decay of the condensate into the contin- 
uum at higher densities makes the fit less reliable and mean- 
ingful, and we have therefore limited the analysis using this 
technique to densities below 10 12 crrT 3 . 



B. Modulation frequency 
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FIG. 5: (Color online) Resonance curve of conversion effi- 
ciency vs modulation frequency for different pulse durations, for 
a thermal 85 Rb gas. Here B av = 156.41 G, E™/h = -5.39kHz, 
Bmod = 0.065 G, T = 20 nK and n = 10" cnT 3 . Inset: Conversion ef- 
ficiency after a 6 ms pulse with T = 20 nK and T = 80 nK, showing 
its weaker dependence on modulation frequency at higher temper- 
atures. The dashed curves are thermal averages of the perturbative 
estimate of Eq. for 6 ms pulses. 

Resonant behaviour was observed in Ref. fTrill in the strong 
variation of the conversion efficiency with modulation fre- 
quency, which is reproduced by our calculations. The con- 
version efficiency from a thermal gas due to a pulse of fixed 
duration and varying frequency is shown in Fig. [5] For the res- 
onance curve representing 6 ms pulses, the full-width at half- 
maximum is 0.75 kHz. From the Lorentzian fit to the 6 ms 
pulse in Fig. 1 of Ref. [J_6| we extract 0.9 kHz. At longer times 
many-body effects may lead to the production of molecules 
by thermalisation. This could lead, for example, to the pro- 
duction of molecules for modulation frequencies smaller than 
—EL^/h, and so increase the width of the resonance curve. The 
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inset of Fig. [5] shows estimates of the conversion efficiency 
using the perturbative estimate of the transition amplitude in 
Eq. (0. The agreement with the numerical result is signif- 
icantly better than that of the transition probability density, 
shown in Fig. [2] due to the effect of thermally averaging over 
all of the continuum states. 

The maximum of a thermal distribution is at a higher en- 
ergy in a warmer gas, and so the optimal modulation fre- 
quency increases with temperature. However, the dependence 
of the conversion efficiency on modulation frequency weak- 
ens at higher temperatures, as shown in the inset of Fig. [5] 
This is caused by the changes in the thermal distribution of 
the gas, which has a decreasing maximum and an increasing 
width as the temperature rises. The decreasing maximum of 
the distribution leads to less being gained by optimising the 
modulation frequency w mo( j/27r. Conversely, the increasing 
width means that a wider range of w mo d have a significant 
population of atoms close to the resonant continuum energy 
Pres/m = E^ v + fi<w mo( j. In general, the stronger resonant be- 
haviour in colder gases allows more efficient conversion. 



5 




999.5 1000 1000.5 1001 1001.5 1002 1002.5 1003 
Modulation frequency (kHz) 

FIG. 6: (Color online) Conversion efficiency vs modulation fre- 
quency for a thermal 133 Cs gas, from a 6ms pulse. The gas 
density and modulation amplitude are identical to the 85 Rb curve 
shown in Fig.[5] Here T = 20nK, n = 1 X 10" cm 4 , S av = 21.37 G, 
B mod = 0.065 G, and E™/h = -1 MHz. 

We have studied the conversion efficiency vs frequency for 
molecular binding energies of h x 100 kHz, and found that it 
leads to a resonance curve of similar width and maximum to 
that for the binding energies examined above. 85 Rb2, though, 
is unstable with respect to inelastic spin relaxation lfl4l l27ll . 
We have therefore also performed the calculation for 133 Cs 
atoms in the (F = 3, rap = 3) Zeeman ground state for a 
molecular bound state energy of —h x 1 MHz. A similar res- 
onance curve is obtained, as shown in Fig. [6] It is broader 
and has a lower maximum than the comparable calculations 
for 6 ms pulses in 85 Rb, which had an identical modulation 
amplitude and gas density. Despite the deeper binding en- 
ergy, the conversion efficiency grows at a similar rate. The 
evolution of the transition probability density for a continuum 



state depends primarily upon its detuning from the resonant 
continuum energy. Consequently, it is primarily the width of 
the thermal distribution, rather than the molecular bound state 
energy, that determines the order of magnitude of the pulse 
duration necessary for association. Our calculations indicate 
that resonant association can be efficient for binding energies 
ranging from h x 5 kHz to h x 1 MHz. It is necessary, however, 
that the chosen molecular binding energy be sensitive to varia- 
tions in the magnetic field. If this is not the case, the magnetic 
field modulation has little or no effect on the diatomic level 
spectrum, and so significant transitions between the contin- 
uum states and the molecular bound state do not occur. Such 
a weak dependence on the magnetic field can occur due to 
an avoided crossing with another bound state, as occurs for 
133 Cs2 at some binding energies J36D . In some cases, it may 
be possible to compensate for this by using a larger modula- 
tion amplitude. 



C. Modulation amplitude 

In the experiments of Thompson et al, increasing the mod- 
ulation amplitude B mo d with a fixed frequency and pulse du- 
ration gave a point of maximum conversion, and after reach- 
ing a minimum a partial revival was observed 13711 . Examin- 
ing the transition probability density of Eq. © for different 
Bmoi shows that as B mo( j is increased, the resonant growth 
of p(p,es,0 is at first amplified, as shown in Fig. [7^. The 
faster resonant growth is also reflected in the proportional- 
ity of the analytic estimate of T(p, f) in Eq. ((T) to B rao d. For 
Bmod = 0.35 G, there is no longer resonant growth in p(p le s, t) 
over a 1 ms pulse duration, although this modulation ampli- 
tude does maximize p{p K>i , 1 ms). The changing amplitude 
and position of the maximum as B mo d varies alters the qual- 
ity of the fit to the thermal distribution, and consequently the 
conversion efficiency. 

For B mo d < 0.9 G, resonant growth is still observed; how- 
ever, the continuum energy of the resonantly growing state is 
negatively shifted from that predicted by Eq. (|5). As shown 
in Fig. |7j3, p(p, 1 ms) has a peak in momentum which, as B mo d 
is increased, at first grows in amplitude and retains its width 
and position, before being shifted towards p = 0. Fully de- 
structive interference for continuum energies up to a few kHz 
occurs when B mo d ~ 1 .0 G and so a minimum in conversion 
efficiency is produced, as reflected in Fig. [8] Beyond this 
modulation amplitude, there is no continuum energy for which 
quadratic growth of p(p, t) is observed. Figure |7p also shows 
two bands of constructive interference in p(p, f). These bands 
have a peak energy which is also dependent on B mo d- Conse- 
quently, at the modulation amplitudes for which these peaks 
coincide with the thermal distribution, a revival of the conver- 
sion efficiency occurs. 

The maximum, minimum and revival in conversion effi- 
ciency are shown for three different modulation frequencies in 
Fig. [8] The absolute conversion efficiency at the maximum is 
strongly temperature dependent, as shown for B mo d = 0.065 G 
in Fig.|5j however, the modulation amplitude giving the maxi- 
mum conversion has a weak variation with temperature above 
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FIG. 7: (Color online) The transition probability density p(p, t) in 
a 85 Rb gas for different values of B moi . Here B av = 156.352 G, 
a> mod /2n = 6.5 kHz, and E^/h = -4.88 kHz. (a) The evolution of 
p(p, t) for the continuum state of energy satisfying the resonance con- 
dition of Eq. ^9). Quadratic growth ceases to be observed when the 
modulation amplitude becomes too great, (b) The transition prob- 
ability density distribution p(p, 1 ms) for B mo *i = 0.15, 0.3, 0.8 and 
1.5 G, as indicated in the contour plot shown in the inset. For each 
fimod < 0.9 G, the peak in energy of p(p, t) grows resonantly on a 
timescale of 1 ms. The continuum energy of the maximal p(p, 1 ms) 
is negatively shifted from p^ cs /m with increasing B rao( i. Two bands of 
revival in p(p, t) can be seen in the inset, as well as in the distributions 
of p(p, t) for S mod = 0.8 and 1.5 G. 



20 nK. The continuum energy for which the state is reso- 
nantly coupled varies with the modulation amplitude, and so 
both frequency and amplitude should be matched to the tem- 
perature of the gas for maximum conversion. Of the plots 
shown in Fig. [8] for example, the best conversion is achieved 
for B mo d = 0.7 G and w mo d/27r = 6.5 kHz. The revival in con- 
version efficiency for larger modulation amplitudes occurs in 
a peak that is narrower, and is due to the constructive interfer- 
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FIG. 8: (Color online) Conversion from a thermal 85 Rb gas as a func- 
tion of modulation amplitude for 6 av = 156.352 G, n = 10 11 cm" 3 , 
T = 20 nK and (±> mo a/27t = 4.9, 5.5 and 6.5 kHz. The solid line 
(<^moi/2n = 6.5 kHz) is a thermal average of the data shown in Fig.|7] 
The variation in conversion efficiency with B rao d is caused by the 
changes in the distribution of p(p, t) shown in Fig.|7p. The revivals 
are caused by the regions of constructive interference, shown in 
Fig Ep, coinciding with the thermal distribution. 



ence shown in Fig.|7J). We note that for B mo d > 1.352 G the 
resonance position Bo is being crossed during the pulse. 



V. CONCLUSIONS 

Resonant association has been experimentally shown to be 
an effective technique of producing molecules llm [l7ll . Here 
we have studied the dependence of the conversion efficiency 
on the duration, frequency and amplitude of the pulse, and 
the density and temperature of the gas. We have shown that 
for a homogeneous gas, the continuum shapes the dynam- 
ics of the association in such a way that it is unlike a two- 
level system, in contrast to the case of resonant association in 
strongly confining optical lattices. The presence of other con- 
tinuum states around that resonantly coupled to the Feshbach 
molecule leads to the requirement of optimising the proper- 
ties of the pulse for the gas in question. Maximum conversion 
requires the amplitude and frequency of the modulation to be 
together optimised for the density and temperature of the gas. 
Colder gases have narrower thermal distributions and so dis- 
play stronger resonant behaviour. The width of the thermal 
distribution also leads to the dephasing of the oscillations in 
conversion efficiency observed at short times in Ref. lfl6Tl . An 
increase in temperature causes a positive shift in the optimal 
frequency for association, but also lowers the maximum pos- 
sible conversion efficiency. 

The amplitude of the modulation and mean-field shifts lead 
to the resonant coupling of continuum states of different en- 
ergy, and thus also affect the conversion efficiency. A higher 
modulation amplitude causes a less energetic continuum state 
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to be resonantly coupled. Beyond a certain amplitude, no res- 
onant growth in transition probability density occurs; how- 
ever, for the parameters of Ref. Il6ll a revival in conversion 
efficiency is seen due to a region of constructive interference 
between the different continuum states. A weak dependence 
of the molecular binding energy on magnetic field limits the 
effectiveness of resonant association, although this can some- 
times be compensated for by an increase in the modulation 
amplitude. The evolution of the transfer probability density 
from a state is primarily determined by its detuning from the 
resonant continuum energy. Consequently, the pulse dura- 
tion necessary for association does not vary significantly with 
the molecular binding energy. We have performed calcula- 



tions for molecular binding energies ranging from h x 5 kHz 
to h x 1 MHz, and predict that resonant association can be ef- 
fective over this range. 
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